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Saving energy and enhancing performance are secular preoccupations shared by both nature and 
human beings. In animal locomotion, flapping flyers or swimmers rely on the flexibility of their wings 
or body to passively increase their efficiency using an appropriate cycle of storing and releasing elastic 
energy. Despite the convergence of many observations pointing out this feature, the underlying 
mechanisms explaining how the elastic nature of the wings is related to propulsive efficiency remain 
unclear. Here we use an experiment with a self-propelled simplified insect model allowing to show 
how wing compliance governs the performance of flapping flyers. Reducing the description of the 
flapping wing to a forced oscillator model, we pinpoint different nonlinear effects that can account 
for the observed behavior — in particular a set of cubic nonlinearities coming from the clamped- 
free beam equation used to model the wing and a quadratic damping term representing the fluid 
drag associated to the fast flapping motion. In contrast to what has been repeatedly suggested in 
the literature, we show that flapping flyers optimize their performance not by especially looking for 
resonance to achieve larger flapping amplitudes with less effort, but by tuning the temporal evolution 
of the wing shape (i.e. the phase dynamics in the oscillator model) to optimize the aerodynamics. 



INTRODUCTION 

Flying animals have since long inspired admiration and 
fueled the imagination of scientists and engineers. Along- 
side biologists studying form and function of flapping fly- 
ers in nature [TJ |2] , the last decade has seen an impressive 
quantity of studies driven by engineering groups using 
new techniques to develop and study artificial biomimctic 
flapping flyers [3, 4 . The widespread availability of high- 
speed video and in particular the merging of experimental 
methods borrowed from fluid mechanics into the toolbox 
of the experimental biologist have permitted to elucidate 
various key mechanisms involved in the complex dynam- 
ics of flapping flight (see e.g. [MI])- 

A recent field of investigation concerns the efficiency of 
flapping flyers, the major interrogation being about how 
natural systems optimize energy saving together with 
performance enhancement. In particular, the passive role 
of wing flexibility to increase flight efficiency through the 
bending of the wings while flapping has attracted a lot of 
attention. It is commonly agreed that this efficiency en- 
hancement comes from the particular shape of the bent 
wing, which leads to a more favorable repartition of the 
aerodynamic forces (see [8] and [9] for an extensive re- 
view). For flying animals in air, such as insects, it has 
been proposed |1UH12| that wing inertia should play a 
major role in competing with the elastic restoring force, 
compared to the fluid loading. The mechanism govern- 
ing the propulsive performance of the flapping flyer can 
therefore be seen at leading order as a two-step process, 
where the instantaneous shape of the wings is determined 
by a structural mechanics problem which then sets the 
moving boundaries for the aerodynamic problem. 

From a dynamical point of view, if we consider chord- 
wise bending of a wing with a given flapping signal im- 



posed at the leading edge, the instantaneous shape of 
the structure is strongly dependent on the phase lag be- 
tween the forcing and the response of the wing (respec- 
tively the leading and trailing edges). Recent works by 
[T3] and [TJ using a simplified model of a flexible wing 
as a combination of heaving and passive pitching have 
shown that a transition from enhanced thrust to under- 
performance occurs for a critical phase value close to the 
resonant frequency of the system. This sustains the com- 
monly invoked argument suggesting that flapping flyers 
could take advantage of a structural property to save en- 
ergy by matching the relaxation frequency of their com- 
pliant wings to the wingbeat frequency |131 rToTUT] . In 
nature this has been observed in particular for undu- 
latory swimming fish or other swimmers that use de- 
forming propulsive structures, such as jellyfish or scal- 
lops (see [TB] and references therein). In the case of in- 
sects, however, the few available observations (especially 
for large species) report wingbeat frequencies far below 
the natural relaxation frequencies [H5H22]- Recent exper- 
iments using a self-propelled model with large-flapping- 
amplitude elastic wings |12j are consistent with the lat- 
ter, since the propulsive efficiency of the model peaks 
for a flapping frequency lower than the primary linear 
resonance of the wings. Fully predicting the wing beat 
rate as the undamped resonant frequency of a linear os- 
cillator (see e.g. [T5]) should be therefore taken with 
reserve. Super-harmonic nonlinear resonances have been 
invoked [23], suggesting that flying animals may effec- 
tively flap their wings far below the primary resonance 
while increasing their performance. This is probably one 
mechanism among others governing the dynamics of flap- 
ping flyers, but it is clear that the details of the underly- 
ing fluid-structure interaction problem are poorly under- 
stood. More specifically, the underlying phase dynamics 
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FIG. 1: Experimental setup: 
Marey |24]. b) Actual setup, 
model used for this study. 
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that set the instantaneous wing shape and lead first to 
an increase and then a loss of the thrust power (and even 
a reversal of the propulsive force as in the case of |13| ) 
remain unexplained. 

In this paper we address these questions using the ex- 
perimental self-propelled flapping-wing model with elas- 
tic wings described in f!2| . Exploring a wide range of 
bending rigidities we show that, in the simplified context 
of chordwisc-compliant wings, the performance optima of 
the system are far from being set by a simple resonant 
condition. We develop a nonlinear one-dimensional beam 
model for the bending wing which is reduced to a forced 
oscillator model suitable to study different nonlinear ef- 
fects. In particular, a set of cubic nonlinearities coming 
from the clamped-free beam equation and a quadratic 
damping term representing the fluid drag associated to 
the fast flapping motion permit to account for the ob- 
served behavior. We show that the nonlinear nature of 
the fluid damping is an essential feature to determine 
the phase lag that leads to an increase/decrease of the 
efficiency. 

As a whole fluid-solid interaction process leading to 
propulsion, we provide evidence that flapping flyers may 
optimize their performance not by especially looking for 
resonance but by using passive deformation to streamline 
the instantaneous shape of the wing with the surrounding 
flow. 



EXPERIMENTS 
Setup and physical quantities 

The experimental setup is the same described by Thiria 
& Godoy-Diana [H], inspired from the pioneer 19 th cen- 
tury experiment by Marey |24| : a flapping wing device is 
attached to a mast that is ball bearing mounted to a cen- 
tral shaft in such a way that the thrust force produced by 
the wings makes the flyer turn around this shaft. A par- 
ticular attention has been paid to reduce friction losses 
in the whole system. Wings are made of Mylar® semi- 
circles of diameter S — 2L = 6 cm. The experimental 
parameters are the forcing frequency (/), the flapping 
amplitude (A u ) and the chordwise rigidity of the wings 
(B) governed by their thickness h. In contrast with the 
first study reported with this setup , the set of wings 
used here covers a larger range of bending rigidities, from 
near-rigid to very soft materials. Six pairs of wings have 
been tested. Their structural properties (thickness, mass, 
and rigidity) are summarized in Table 1. 

This specific setup allows to measure various averaged 
quantities (see [T5] for details): the cruising speed U 
when the device is allowed to turn around, and the thrust 
force Ft when it is held at a fixed station (see Fig. [2] (a) 
and (b)) which gives the averaged aerodynamic thrust 
power, being the product Pt — UF. In both cases, the 
power consumption Pj is measured. On the other hand, 
we performed a precise dynamical study of the flapping 
wing. For each set of parameters (A u , ff, B), the phase 
and amplitude of the trailing edge, with respect to the 
forcing flapping motion, has been measured using a fast 
cadenced camera (1000 fps) in both air and vacuum. 

It is important to recall that for this setup, and more 
generally for flapping flyers in air, the main bending mo- 
tor of the flexible wings is wing inertia |10lfl2] . The com- 
petition between the wing inertia and the elastic restoring 
force is captured by the scaled elasto-inertial number Af e i 
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The first expression is a direct comparison between 
both the moments of inertial and elastic forces. Inter- 
estingly, this number can also be expressed as a func- 
tion of the ratio between the forcing and relaxation fre- 
quencies times the non-dimensional forcing amplitude of 
the driving motion, which allows to express directly the 
bending rate as function of a non-dimensional oscillator 
forcing term. The second expression is therefore useful 
to explore the nearness of the resonance and will be used 
to analyze the experimental data in this paper. Results 
will therefore displayed as a function of the reduced fre- 
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FIG. 2: (a) Cruising speed, (b) thrust force and nondimen- 
sional (c) thrust (pr) and (d) input (pi) powers as a function 
of uJj. The gray area represents the optimum region, the 
dashed line indicates the location of the reduced natural fre- 
quency of the wing (linear resonance). 



the reduced flapping amplitude. In order to compare the 
aerodynamic performance in all the experiments, both 
the thrust force and cruising speed were rendered non- 
dimensional using the appropriate scalings fx = FtL/B 
and u = U/A^ui. The non-dimensional powers (displayed 
in Fig. [2](c) and (d)) then read pr — UFtL/BA^lu and 
Pi = PiL/Buj. 

In both the thrust force and cruising speed curves, it 
is clear that increasing wing flexibility brings out two 
distinct regimes: up to a certain flapping frequency, the 
more flexible wings outperform the rigid linear U(f) re- 
lationship (see also [25 ). The measurements for the two 
most flexible wings evidence the appearance of an under- 
performance regime in which both Ft and U lie below 
the rigid wing case. Looking now at the nondimensional 
thrust power, the data from all wings collapse on a single 
curve with a clear performance peak, which agrees with 
what has been observed by [131 [Hj for heaving/pitching 
systems. An important point is that the maximum in 
performance does not take place at the resonant fre- 
quency, but much below (around 0.7wo, represented by 
the gray shaded area). Moreover, the nondimensional 
thrust power at cjy — 1 (see dahsed line in Fig. [2] (c)) is 
even more than 4 times lower than the optimum value. 
At last, we remark that there is also no sign of a resonant 
behavior in the consumed power curve (Fig. [2] (d)). 



Wing dynamics 

We proceed now to study the behavior of the wings 
considered as a forced oscillator, assuming the oscillation 




(b) 



0.02 0.04 0.06 0.08 













' A. 















0.1 



0.02 

f(s) 



0.04 



FIG. 3: a): Photograph of the flapping wing showing succes- 
sive states of the bending wing during one stroke cycle (thick- 
ness is 0.050 mm and cJf — .5). As can be seen, the main 
deformation is mainly performed on the first mode. In this 
case the phase lag is quite large, leading to a strong increase 
of flight performance, b): Typical time series tracking the 
motion of the leading (black curve) and trailing (red curve) 
edges of the wing at mid-span, obtained from video record- 
ings at 1000 fps. c): Same as b) but with a forcing near |wo, 
exhibiting super-harmonic resonance typical from dynamical 
systems containing cubic nonlinearities. 



of the leading edge to be the forcing and that of the 
trailing edge to be the response (which means to assume 
that the wings bend following only the first deformation 
mode). As said before, the amplitude and phase shift 
of the response can thus be measured by following the 
two wing edges on a high cadenced camera recording (as 
seen on Fig. [3] (a)). Figs. [3] (b) and (c) display two 
characteristic time evolutions of the driving oscillation 
(the imposed wing beat, shown as black dots) and the 
wing elastic response (the motion of the trailing edge, 
red dots) in the moving frame. The first case shows a 
typical response, at ulf = 0.79, mainly sinusoidal at the 
driving frequency, which supports the assumption that 
the oscillations of the wing follow a single mode. In the 
second case, the driving frequency is near one third of 
the resonant frequency ljq. As can be observed in Fig. [3] 
(c), the response is then a combination between ojo/3 and 
wo, giving evidence of a super-harmonic resonance |26| . 
pointing out the fact that the system integrates cubic 
nonlinearities. The non dimensional amplitude a (i.e. 
scaled by the length of the wing L) and phase 7 have 
therefore been extracted from those signals for each pair 
of wings as a function of the reduced driving frequency for 
two different amplitudes. Results are displayed in Fig. [4] 
In parallel, the same experiments have been conducted 
in a vacuum chamber at 10 % of the ambient pressure. 
Results are also displayed in Fig. [4] for comparison. 
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Table 1. Wing properties 



wing thickness, h (mm) 


0.050 


0.078 


0.130 


0.175 


0.250 


0.360 


mass per unit area /i 3 (kg.m -2 ) 


4.50 10" 2 


10.63 10~ 2 


17.67 10" 2 


24.12 10" 2 


34.92 10" 2 


47.95 10" 2 


rigidity B (N.m) 


3.34.10 -5 


1.83. 10 -4 


1.02.10~ 3 


2.26. 10 -3 


7.31. 10 -3 


14.00. 10 -3 


relaxation frequency fo (Hz) 


25.4 


34.2 


62.2 


89.5 


117.1 


160.8 


color label in figures 


blue 


red 


green 


yellow 


purple 


black 



As can be seen, the evolution of the amplitude a shows 
a fast increase from very low flapping frequencies. This 
is the expected behavior owing to the inertial charac- 
ter of the forcing. A slight but rather broad peak can 
be observed in the nearness of u; /3 in the amplitude 
curve, confirming the occurrence of the super-harmonic 
resonance hinted above and strengthening the fact that 
this type of mechanism may play a role as a strategy 
for performance enhancement in nature |23j . Two more 
points have to be underlined: first, measurements in air 
and vacuum are approximately the same, in accordance 
with the hypothesis that inertia is the main bending fac- 
tor for flapping flyers |10H12| . The second point is that 
no clear resonance is observed around ulf — 1 (only a 
barely visible peak in the case of the lowest forcing am- 
plitude shown in the insert in Fig. |4ja)). Concerning the 
phase 7, the present results recover the trend of what 
has been observed recently HS1 UH HS] : M increases 
monotonically with tJf. Considering the experiments in 
air at normal conditions, this observation together with 
the performance increase shown in the first part of the 
Pt(&f) (Fig. [2](c)), brings the following conclusion: the 
increasing phase shift 7, which corresponds to a situation 
where the wing experiences a larger bending at the max- 
imal flapping velocity, leads to a more favorable reparti- 
tion of the aerodynamic forces (as discussed in |12j). 
A simple argument widely shared in the community con- 
necting the phase dynamics to the propulsive perfor- 
mance is: the larger the phase lag is, the best the thrust 
power would be |13[ 114] , until the point where the wing 
experiences its largest bending at 7 — tt/2. However, 
while the argument reasonably agrees with the observa- 
tions in the range of forcing frequencies where perfor- 
mance increases with uif, the maximum performance does 
not actually match with the maximum of bending that 
occurs at 7 = tt/2, but relatively far below this expected 
optimum (which lies actually around ir/4). 
One last important remark to be made concerns the 
phase evolution in vacuum. It is clearly observed that 
7 decreases more slowly in the low density environment 
within the whole range of flapping frequencies studied. 
In contrast with the amplitude measurements, where 
the data from the experiments in vacuum follow roughly 
the same curve of those in air at atmospheric pressure, 
the large difference in the 7 curves between both cases 
points out unequivocally the importance of the surround- 
ing fluid in determining the phase dynamics. This point 
will be discussed later. At this stage, we have shown that, 
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FIG. 4: Evolution of the non-dimensional amplitude a) and 
phase b) of the trailing edge wing response as a function of the 
reduced driving frequency for both flapping amplitudes A u = 
0.8 and A w = 0.5 (filled symbols correspond to measurements 
in air, open symbols in vacuum). T hose results are compared 
to nonlinear predictions from Eq. 12 with (gray line) and 



without (black line) nonlinear air drag (discussed further in 
the text). 



as observed in the pitching/heaving systems of [T31 114j . 
the increase in performance of elastic wings undergoing 
large oscillations is essentially governed by a fast grow- 
ing phase evolution. However, the physical mechanisms 
governing the propulsive performance remain unclear. In 
particular, the mechanisms leading to the useful evolu- 
tion of 7 as well as the link between resonance and per- 
formance are still looking for a definitive answer. 



NONLINEAR ID BEAM MODEL 

In order to understand those crucial points, one can 
consider the elastic wing as a clamped-free beam under 
base harmonic forcing. For simplicity, the beam is con- 
sidered as one-dimensional taken at mid length in the 
spanwise direction of the wing. We assume here, accord- 
ing to the experiment, that only flexural displacements 
(i.e. perpendicular to the direction of the flight motion) 
are allowed. The structural properties of the beam are 
determined by measuring experimentally the relaxation 
frequency. 

Thus, the equation governing the motion of the non- 
linear flexural oscillations of clamped-free beam writes 
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where is the transversal local displacement, the 
Young modulus. / the second moment of inertia and fi 
the mass per unit of length. Writing W as W(x,t) = 
w(x, t) +wo(t), where wo(t) is the driving motion defined 
by Wo(t) — Aucosfaft), and using the non-dimensional 
quantities for space and time w 



t = (^f) 1 2 L 2 , equation |2j reads: 
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which has to satisfy the clamped-free boundary condi- 
tions w(0,t) = w'(0,i) = w"{l,t) = w"'{lj) = 0. The 
last term on the right hand side in Eq. pi — A u w$ = 
A^olf 2 cos(u7/t) = J\f e i cos(cJ/f), is a forcing term due to 
the wing inertia whose amplitude is given by the elasto- 
inertial number and which is dependent on the square of 
the driving frequency as seen before. 
The next step is to set apart the spatial dependence by 
projection of Eq. [3] onto the complete set of eigenfunc- 
tons defined by the linear part. The displacement is ex- 
pended as w(x,t) — X p (t)$> p (x) (see [28J) where <j> p 
are the non-dimensional linear modes for clamped-free 
beams which are not recalled here for the sake of brevity. 
The problem then writes (the ~ have been removed for 
simplicity) : 
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where h p - k and ffj k are determined by: 
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The projection of the forcing term on the p th mode. 
F pi writes at the trailing edge: 



F p = A u wf%(l) / $p(x)dx 



(7) 



As the propulsive regimes observed in this work lie be- 
low the first relaxation frequency of the wing, we assume 
that the response of the wing is mainly governed by the 
first eigcnmode. Hence, equation [3] can be considerably 
simplified and reduces for the only mode 1 to: 



X + X = -h\ n X 3 — fi U (X 2 X + XX 2 ) + F 1 (t) (8) 

A crucial feature is now to choose a damping term 
to this dynamical system. During a stroke cycle, the 
wing follows very fast motions involving high local 
Reynolds numbers, which prompt us to include a non- 
linear quadratic fluid drag term [29 in addition to the 
classical linear viscous friction law. The damping is then 
chosen as a combination of linear and nonlinear terms as 
follows: 



E(X,X)=SX + Z lU \X\X 



(9) 



The linear and nonlinear coefficients £ and £ n ; are 
estimated studying the impulse response for each wing 
(21)] The solution of Eq. [8] including damping is deter- 
mined by using a classical multiple scale method at first 
order (see [26 ). To this end, we introduce a small pa- 
rameter e and a detuning parameter a — (oJy — l)/e. The 
problem to be solved reads. 



X + X = -e(h{ n X 3 + fl n (X 2 X 
+ ~(X,X)+F 1 (t)) 



XX 2 ) 



(10) 



According to the multiple scales theory, we express 
the solution in terms of different time scales as X — 
X (to,ti) + eXi(t ,ti) + .... where t — t and t\ = et 
are respectively short (relative to the oscillation of the 
wing) and long times scales. The system at order e° is 
d 2 Xo + Xq = an gives the straightforward solution 
X = A(ti)e ita + A*(t 1 )e~ ita where A and A* are com- 
plex functions. 

At order e , we obtain: 



X\ — — h\ n XQ — / 1 1 11 (AqA + X Xq ) 

- E(X , X ) - 2d tlto X + F x cos(i + trti) 

(11) 



Using the expression of X found at order e° into Eq. 
1 1 1 1 an equation for A is obtained by elimination of the 
secular terms: 
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A 2 A*(3h\ 11 ~2fl 11 )+z(2d tl AHA+^\A\A) = l -F ie 

(12) 

where the pre-factor ~ in front of the nonlinear damp- 
ing coefficient is obtained during the special integration 
over one period of the Fourier expansion of the function 
Xo|Xq| (see [25]). As can be seen. Eq 



12 



is a charac- 
teristic equation of a forced damped oscillator with cu- 
bic nonlinearities. At last, substituting the polar form 
A = iae 1 ^* 1-7 ) , separating into real and imaginary parts 
and looking only to the steady-state solutions, we find 
two relations for the amplitude a and phase 7. [32 





A„/L = 1.2. 




2^2 



Ff 



7 = arctan 



Via 3 — aa 



(13) 
(14) 



where T\ = ^(3h\ 11 — 2 flu) is the nonlinear cubic term 
coefficient, which is computed from Eq. [5] and [6] 

Eq. [14] closely resembles a classic nonlinear Duffing 
oscillator except that the forcing amplitude is frequency 
dependent and that a nonlinear damping term is present. 



ANALYSIS AND DISCUSSION 



Resonance and phase evolution 



FIG. 5: Dependence of the amplitude a and phase 7 with 
the reduced forcing frequency uJf for the first mode of a 
clamped-free beam forced by inertia for two different (high 
and medium) amplitudes A u (chosen arbitrarily for clarity). 
The blue line corresponds to the linear prediction, the black 
line to the non-linear model from Eq [12] with linear damping 
— 0), the red line to the non-linear model with nonlinear 
damping (£ n ; 7^ 0). As can be seen only cases with relatively 
small flapping amplitude and linear damping can exhibit a 
slight resonance peak. Greater amplitudes and/or presence 
of nonlinear damping behave as a non-resonant system in the 
domain of flyers capabilities. Concerning the phase, models 
including only linear damping do not produce "useful" phase 
lag except in the nearness of the phase jump. In contrast, the 
presence of a nonlinear damping produces a fast and helpful 
evolution. 



Predictions of the above model for the parameters of 
the experiments are plotted in Fig. [4] for both cases in 
air and vacuum. In addition, for a clear understanding 
of the underlying dynamics described by Eqs. [T3| and [T4] 
a comparison between predictions from a linear model, 
a nonlinear with linear damping and a nonlinear with 
nonlinear damping is displayed in Fig. [5] for two flapping 
amplitudes . It can be seen that the model based on a 
single mode is capable of reproducing all the observations 
made from the experiments both in normal and low den- 
sity environments. The good agreement between experi- 
ments and model allows us to pinpoint some mechanisms 
underlying the complex mechanisms of flapping flight. 
The first concerns the question of resonance: from Fig. [5] 
it can be observed that the only case (apart from the lin- 
ear case) exhibiting a slight resonance peak corresponds 
to relatively small flapping amplitude and damping co- 
efficient [i.e. only linear damping term, see Fig. [5] (a)]. 
Cases for higher amplitude and/or presence of nonlin- 
ear damping behave as a non-resonant like system in the 
range of flapping frequencies studied. In nonlinear oscil- 
lators, it is known that the main effect of the nonlinear 
term is to distort the resonance curve and shift the reso- 
nance peak to higher frequencies (for a hardening coeffi- 
cient Ti > 0, as in the present study) [26 . An important 



feature of such nonlinear systems is that the distortion of 
the shape of the resonance curve is directly dependent on 
the amplitude of the excitation. In the present case where 
the forcing is inertial, the response depends on the square 
of the forcing frequency (or on the elasto-inertial number 
M e i), which provides an increase of the amplitude plotted 
in Fig. [4] independent of an intrinsic resonance mecha- 
nism. Hence, we can expect the actual resonance curve 
of the system to be all the more distorted that the flap- 
ping frequency increases. Another feature that makes it 
difficult for the flapping flyer to benefit from a resonance 
mechanism is the presence of a geometric saturation due 
to the finite length of the wing. Always due to the in- 
ertia effects, this geometrical saturation will be reached 
all the more soon that the demand for larger amplitude 
(i.e. better performances) is increased. Coming back to 
the distorted resonance curve, the visible consequence is 
that the wing, even for a small nonlinear cubic coefficient, 
behaves as a system never reaching a peak in the range 
of frequencies commonly used by flapping flyers. Addi- 
tionally, the presence of strong damping accentuates this 
behavior by smoothing the value of a possible resonance 
peak. This last observation is consistent with the fact 
that birds or insects may not especially look for struc- 
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tural resonance to improve their performance. 

The second point is the crucial role of fluid damping 
in triggering the phase lag that is useful for thrust en- 
hancement. For the phase, shifting the resonance peak 
as a result of the nonlinear spring in the oscillator model 
means shifting the phase jump at 7 = tt/2 to higher 
frequencies as well. Thus, without air drag, as can be 
seen in Fig. [5] (c) and (d), the nonlinear evolution of 
the phase 7(0)/) would be even slower than in the linear 
case for which the phase evolution is already not espe- 
cially favorable except in the nearness of the resonance. 
This is exactly what is observed for the vacuum measure- 
ments where the nonlinear damping due to fluid drag is 
negligible. On the contrary, the presence of a quadratic 
fluid damping determines a fast increase of the phase lag 
(and a so a thrust improvement) even from the very first 
flapping frequencies. This implies of course that strong 
flapping velocities are a necessary condition for the bend- 
ing to become efficient (i.e. elasticity will play a minor 
role if the flapping beat amplitude is not strong enough) . 
Summarizing, the instantaneous wing shape is given by 
the two following ingredients: inertia provokes the bend- 
ing (gives the amplitude) and damping, by controlling the 
phase lag, allows this bending to be usefully exploited. 
Large phase lags will provide largest bending of the wing 
at maximum flapping speed, leading to a more favorable 
repartition of aerodynamics forces. 

Optimum 

Since classic resonance mechanisms cannot answer it, 
the question of the performance optimum (or the tran- 
sition to underperformance) remains unclear. We there- 
fore proceeded to study the kinematics of the wing in 
the laboratory frame. In particular, we have compared 
both characteristic angles relative to the global wing 
motion. The first characteristic angle is dependent on 
the ratio between the maximal vertical flapping veloc- 
ity u u — ujAu and the cruising velocity U and reads: 
<f> = arctan(wA w /t/). is considered as the instanta- 
neous angle of attack of the wing and as can be seen, 
is directly related to the Strouhal number St = ujA^/U 
which determines as well the performance of flapping fly- 
ers |30| . We define a second characteristic angle 8 as 
the trailing-edge angle taken at the maximum flapping 
velocity. This angle is directly related to the phase lag 
7, and thus determines to what extent the bending of 
the wing will be useful in terms of performance. Fig. [6] 
shows the evolution of the ratio 8 '/ 1 which is naturally a 
growing function of oJ/ because both an increase in 8 or 
a decrease in <j) lead to an enhancement of the propulsive 
performance. 

The interesting point is that the location of the per- 
formances/under performances transition takes place at 
8/(f> = 1 (i.e. when both angles point instantaneously at 
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FIG. 6: Evolution of the two characteristic angles of the wing 
motion 6 and as a function of the reduced driving frequency 
uff. Two regimes can be distinguished: (I): cj> < correspond- 
ing to the performances increasing stage due to a useful phase 
lag. (II): (f> > 9 corresponding to the transition to under- 
performances due to a loss of the effective wing area. The 
optimum occurs therefore when <f> and point at the same 
direction (best phase lag). 



the same direction). Thus, the optimum value of 8 does 
not corresponds to the maximum bending experienced by 
the wing (which would be the optimal solution) but to 
the moment when the deflection angle matches the angle 
of attack as sketched in Fig. [6] For a rigid wing, because 
8 is fixed (= 0), the optimization problem is here nonex- 
istent and thrust only depends on the driving frequency 
(for a given amplitude). With flexibility and according 
to what has been previously observed, 8 starts increasing 
and tends to align the wing trailing-edge with the flow. 
As discussed earlier, this leads to a more favorable repar- 
tition of the aerodynamics forces as sketched in Fig. [6j 
However, this argument is only valid if the surrounding 
flow is totally attached to the wing (i.e. separation oc- 
curs only at the trailing-edge). A situation where 8 > <fi 
is strongly subjected to flow separation before the wing 
trailing edge. In this case the effective surface relative to 
the aerodynamic load can be expected to be drastically 
reduced leading to a loss of aerodynamic performance. It 
has to be noticed that the value of tt/2, or more gener- 
ally values of phases greater than 8 opt observed in this 
experiments should be, theoretically, more optimal (i.e. 
should give more optimal bending shapes for useful pro- 
jection of forces). However, if a separation occurs, the 
corresponding loss of thrust force (and so cruising speed) 
will accelerate the decoherence of both angles and hence, 
will provoke the subsidence of the performance, as has 
been observed on Fig. [2] The more economic strategy 
to fly is therefore to set 8 <f> which corresponds to the 
optimum way to transfer useful momentum. 
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CONCLUDING REMARKS 

In this work, we aimed at describing the dynamics gov- 
erning the performance of flapping flyers. Considering 
large flapping amplitude and relatively large wings (as for 
big insect species), we have shown that nonlinear and in- 
ertia effects, together with geometric limitation, question 
the prevailing idea that energy-saving strategies in flap- 
ping flight must be related to resonance mechanisms. In 
search of improving performances, animals may actually 
stay below the resonance point. Besides, the nonlinear 
nature of air drag (which implies sufficiently strong flap- 
ping amplitudes) seems to be a fundamental ingredient 
to create the phase lag between the leading and trailing 
edges of the flapping wing that allows the elasticity en- 
ergy to be used at its best. One last comment is that 
the presence of structure resonances for flyers in nature 
is not invalidated by the mechanism described here. For 
instance, small insects may not use much elasticity and 
bending because either their wings are too small or the 
local Reynolds number is not sufficiently high to produce 
enough damping, and thus a useful phase lag. However, 
studies containing a large bank of comparative resonant 
frequencies and wingbeats of insects or birds being rare in 
the literature, it is consequently hard to draw any conclu- 
sion about the existence of two distinct strategies at this 
state. According to biologists, resonant mechanisms lie 
at the muscle level more than in the wing structure itself 
(see [II [31] and reference therein) which would strengthen 
that there is no reason, a priori, for flapping flyers to look 
for structural resonance of the wing. Further analysis on 
such a way would certainly help to discern if there are, 
or not, universal characteristics for flapping flyers. 
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Tardy for her careful reading of the manuscript. This 
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